Hole Structures in Nonlocally Coupled Noisy Phase Oscillators 
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We demonstrate that a system of nonlocally coupled noisy phase oscillators can collectively exhibit 
a hole structure, which manifests itself in the spatial phase distribution of the oscillators. The phase 
model is described by a nonlinear Fokker-Planck equation, which can be reduced to the complex 
Ginzburg-Landau equation near the Hopf bifurcation point of the uniform solution. By numerical 
simulations, we show that the hole structure clearly appears in the space-dependent order parameter, 
which corresponds to the Nozaki-Bekki hole solution of the complex Ginzburg-Landau equation. 
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One representative class of coupled oscillator systems 
is the coupled p/iase oscillators [J 0, S, 0] • For example, 
globally coupled phase oscillators, such as the Kuramoto 
model, have attracted the attention of many researchers 
for a long time d, H, 0] ■ Recently, while coupled phase 
oscillators on complex networks have been investigated 
widely [1, [1, [l^, [ll[ , nonlocally coupled phase oscillators 
have also been studied, which exhibit a remarkable class 
of patterns called chimera where phase-locked oscillators 
coexist with drifting ones [12, Jji, JJ,, X5, 16]. 

In this paper, we demonstrate that nonlocally cou- 
pled noisy phase oscillators can collectively exhibit a hole 
structure in their phase distribution. After briefly re- 
viewing several results on the phase model [H, [l3] , we 
will present our new findings obtained from numerical 
simulations of the Langevin-type equation and its cor- 
responding nonlinear Fokker-Planck equation describing 
the phase model. We will show that a hole structure 
clearly appears in a properly defined order parameter un- 
der suitable conditions, and compare it with the Nozaki- 
Bekki hole solution of the complex Ginzburg-Landau 
equation, for which the modulus displays regions of local 
depression. 

A system of nonlocally coupled noisy phase oscilla- 
tors is described by the following Langevin-type equation 
(LE) for a phase (j){x, t) at location x and time t: 

/OG 
dx' G{x~ x') T{(t){x, t)~(l>{x',t)) +^{x, t). 
-OQ 

(1) 

Here the first term w represents the natural frequency 
common to all the oscillators, the second term the non- 
local coupling, and the last term the additive noise. The 
phase coupling function r(i^), which is a 27r-periodic 
function of 0, satisfies the in-phase condition, i.e., 
dr{(t))/d(t)\^=o < 2]. The spatial coupling function 
G{x) is given by 



G(x) = ^exp(- |a;|) , 



(2) 



which is normalized in the infinite domain. The noise is 
assumed to be Gaussian- white, whose statistics are spec- 
ified by 

(x, t)) = 0, {x, t) C {x', t')) = 2DS {x - x') S{t- t') . 

(3) 

Equation ^ can be derived from a certain class of 
reaction-diffusion systems under suitable conditions, us- 
ing the phase reduction method after adiabatically elim- 
inating a highly diffusive chemical component [15]. 

In Refs. tl5|,|17[, it was shown that the LE ([T]) is equiv- 
alent to a single-oscillator nonlinear Fokker-Planck equa- 
tion (FPE) given by 



dfitp,x,t) 
dt 



d_ 
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V{i>,x,t)f{^,x,t) 



-D 



d^fiip,x,t) 



(4) 



V{ij,x,t)=uj+ / dx' G{x- x') 



X / d^'r{,p-ip')fiij',x',t), (5) 
/o 



where f{ip,x,t) is a space-time-dependent single phase 
distribution function ofip, i-e., the normalized probability 
density that 4'{x,t) takes a value tJj (see also Refs. p^. 

The phase model ([T]) is capable of sustaining traveling 
waves below a critical noise strength, D = Dc, where 
the uniform solution of the FPE ^ undergoes a Hopf 
bifurcation [l^, [l3| . In the vicinity of this Hopf bifur- 
cation point, we can derive a complex Ginzburg-Landau 
equation (CGLE) 

dtA{x, t) ^ {Dc -D)A + ddlA - g \A\^ A, (6) 

from the FPE (ffl by applying the center-manifold reduc- 
tion method [2|, where we introduced the complex am- 
plitude A{x,t) representing the fluctuation of f{ip,x,t) 
in the phase direction as 
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f{ip,x,t) = — + —(A{x,t)e'^'f'+^^^*+A*{x,t)e-'^'''- 
27r 27r V 



iSlt 



(7) 
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Here A is the wavenumber of the phase fluctuation, 
the Hopf frequency, and A* the complex conjugation of 
A. The parameters of the CGLE ([6]) are given by 



Dr. 



9 = 



(8) 



(9) 



where A = arg max/ QTi/l, and F/ is the Fourier compo- 
nent of the phase coupling function defined by 



neighboring oscillators, i.e., L = N Ax, and the periodic 
boundary condition is imposed. In the continuous limit, 
N ^ oo with L fixed, the Langevin simulation would 
exactly correspond to the Fokker-Planck simulation, but 
we can use only finite values of N in actual Langevin sim- 
ulations. Thus, we must consider the effect of finite-size 
fiuctuations, which comes from the finiteness of the oscil- 
lator number within the nonlocal coupling range. Apply- 
ing the finite-size scaling argument developed in Ref. [26| 
(see also Ref. [13]) to the space-time- averaged modulus 
{R) of the order parameter, we can obtain a scaling form 



F(V')- £ Tie 



(10) 



In what follows, we restrict ourselves to the case of A = 1 
and positive Sftg, i.e., the case that the first Fourier com- 
ponent of the phase distribution function has the largest 
critical noise strength, and the Hopf bifurcation of the 
uniform solution is supercritical. 

Since we assume A = 1, Eq. ([5]) and Eq. ^ depend 
only on the first and the second harmonics of the phase 
coupling function, F-i-i and In this case, without 
loss of generality, the phase coupling function can be ex- 
pressed as 



Y{ip) = - sin {ilj + a) +u sin {2ip + f3) . 



(11) 



Now let us introduce a space-time-dependent complex 
order parameter with modulus R{x,t) and phase Q{x,t) 
as 

/>OC 

R {x, t) e*®(^'*) = / dx'G{x- x') e^-^l^''*) 



2tt 



dx'G{x-x') I d^' e''''' f{ilj',x',t) 
1 Jo 



dx'G{x-x')A* {x\t)e-'^* 
~ A*(a;,i)e-'"*. (12) 

In deriving the last expression, we utilized the fact that 
the spatial characteristic length of the complex ampli- 
tude becomes sufficiently long compared to the nonlocal 
coupling length near the critical point [l^, Thus, 
the order parameter corresponds to the complex conjuga- 
tion of the complex amplitude, which is governed by the 
CGLE ©. As is well known, CGLE admits the Nozaki- 
Bekki hole solutions ^23, [HI, [Si, [H [H HHj , for which the 
existence of the amplitude degree of freedom is crucial. 
In the following, we will show that the phase model ([1]) 
can exhibit a hole structure in its space-time phase distri- 
bution, despite its lack of apparent amplitude variables. 

First of all, let us identify the Hopf bifurcation point 
predicted from the FPE (|4]) by numerical simulations of 
the LE ([T]). In the numerical simulations, our continuous 
media of size L is replaced with a long array of N oscil- 
lators with sufficiently small separation Ax between the 



(13) 



where is a scaling function depending on N and {Dc — 
D) only through the combination N^^^{Dc — D). We fix 
our system size L to be 8.0 and vary the number N of 
the oscillators, i.e., L = N Ax — 8.0. For this Langevin 
simulation, we fix the parameter values as 



a = 0.5, u = 0.0, 



which yield the critical noise strength 



(14) 



L'c = cos (0.5) /2, (15) 

and the two essential parameters of the CGLE [^.[23}. 

ci = 3d/5Rd = tan (0.5) , ca = ^g/^g = - tan (0.5) /2. 

(16) 

With these parameter values, spatially uniform oscilla- 
tions are realized in the Langevin simulation, as expected 
from the phase diagram of the CGLE [25|. Figure[T]sum- 
marizes the numerical results in the rescaled variables, 
where N^''^{R) is plotted as a function of N^/'^{Dc - D) 
for several values of N . All curves collapse onto a sin- 
gle identical curve after rescaling, which agrees well with 
the prediction based on the FPE (|4]). This gives good 
evidence for the existence of the Hopf bifurcation at 
D = Dr. 

Now we demonstrate that the phase model (jlj can ex- 
hibit a hole structure. We first carry out a Fokker-Planck 
simulation. The parameter values are chosen as 

a = arctan (0.5) , /3 = 2.3, m = 1.8, (17) 

which give 



Dc = cos (arctan (0.5)) /2, 



(18) 



ci EE ^d/m = 0.5, C2 EE ~ 2.0. (19) 

The corresponding CGLE (O has stable hole solutions 
with these parameter values ^25j . In our numerical simu- 
lations of the FPE ^ , periodic boundary conditions are 
imposed for both phase and space x. The intervals of 
the phase and the space are 27r and L = 102.4, respec- 
tively. The pseudo-spectral method with M = 32 modes 
is applied for the phase. The number of the spatial grid 
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d6/dx = s tanh [b {x — c)] 



(23) 



N"^ (Dc - D) 



FIG. 1: Dependence of the order-parameter modulus on the 
noise strength, obtained from the Langevin-type equation (HI 
using different values of the number of oscillators A'^. The 
data are plotted using rescaled variables, revealing the scaling 
function given by Eq. (|13|l . 



points is = 512, i.e.. Ax = 0.2. Our numerical re- 
sults are unchanged if we further increase the number of 
modes M, the number of grid points TV, or the system 
size L. The initial condition is given by 

f{ip,x,t = 0) = C^[2 + B{x) exp(iV) + S*(a;) exp(-#)], 

(20) 

where Cx is a normalization constant, and 



Bix) 



A{x-L/A)/L 
exp[27ri(a;-i/2)/L] 



(0 < a; < L/2), 
{L/2 <x <L). 



(21) 



Note that this B{x) has one "phase singularity" and sat- 
isfies the periodic boundary condition. The noise inten- 
sity is chosen as D/Dc = 0.9. 

Figure [2] displays the spatial profile of the phase dis- 
tribution function f{il},x,t) obtained from a numerical 
simulation of the FPE (jH). Figure [3] displays the spa- 
tial profile of the order-parameter modulus (a) and the 
phase portrait of the order parameter (b), which are ob- 
tained from Fig. [2] using the relation given in Eq. (H^). 
We can confirm that a hole structure actually appears 
in the order parameter. The small bump at the right 
of the hole is a shock structure due to the collision of 
counter-propagating plane waves emitted from the hole 
structure. Once such a hole structure is formed, it sta- 
bly persists throughout our numerical simulation. This 
hole structure can be well fitted by a non-propagating 
Nozaki-Bekki hole solution ^23^ in the form 

W-a {x) = a tanh [b{x ~ c)] exp [ie{x)] , (22) 



where a, h, c, and s are real parameters. These pa- 
rameters are estimated as app — 0.45, 6fp — 0.23, and 
spp — 0.15 from the numerical simulation of the FPE ([U, 
while the theoretical values for the reduced CGLE ([6]) are 
given by ogl ^ 0.58, 6gl ^ 0.27, and sql ^ 0.12. The 
agreement between the simulation and the theory seems 
reasonable, in consideration of various approximations 
used in deriving the CGLE §^ from the FPE g]) [29]. 

We also carried out Langevin simulations with peri- 
odic boundary conditions. We prepared an appropriate 
initial distribution of the oscillators using the hole solu- 
tion of the FPE (ID) . Figure 2] displays the snapshot of 
the local oscillator phase obtained after the initial tran- 
sient for L = 102.4 and N = 2^"^. The spatial profile of 
the local oscillator phase well corresponds to the spatial 
phase distribution shown in Fig. (5] However, this phase 
distribution corresponding to the hole structure eventu- 
ally collapses to that corresponding to a plane wave due 
to the finite-size fluctuation that we mentioned above. 
The lifetime of the hole structure clearly increases with 
N, so that the hole structure is expected to exist stably 
in the iV ^ oo hmit [soj. 

In summary, we studied a system of nonlocally coupled 
noisy phase oscillators based on the LE ^ and its cor- 
responding FPE ([4]) . We confirmed that the onset of the 
coherence in the order parameter of the LE ^ is identical 
to the Hopf bifurcation of the FPE using the finite- 
size scaling relation for the numerical data obtained from 
the Langevin simulation. We then demonstrated that a 
stable hole structure can appear in the order parame- 
ter calculated from the FPE ([l]), which corresponds to 
the non-propagating Nozaki-Bekki hole solution of the 
CGLE ([6]). Phase models generally lack the amplitude 
variables, which are crucial for the hole solutions with 
phase singularities. In the noisy phase model ([T|), how- 
ever, the external random force effectively produces the 
amplitude degrees of freedom in the phase distribution of 
the oscillators. The phase model ([T]) can exhibit yet an- 
other interesting phenomenon called noise-induced turbu- 
lence [15, 17]. Detailed numerical and theoretical analysis 
on this phenomenon is recently reported in Ref (28| . 
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